Resummation of the Divergent Perturbation Series 
for a Hydrogen Atom in an Electric Field 
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We consider the resummation of the perturbation series describing the energy displacement of a 
hydrogenic bound state in an electric field (known as the Stark effect or the LoSurdo-Stark effect), 
which constitutes a divergent formal power series in the electric field strength. The perturbation 
series exhibits a rich singularity structure in the Borel plane. Resummation methods are presented 
which appear to lead to consistent results even in problematic cases where isolated singularities or 
branch cuts are present on the positive and negative real axis in the Borel plane. Two resummation 
prescriptions are compared: (i) a variant of the Borel-Pade resummation method, with an additional 
improvement due to utilization of the leading renormalon poles (for a comprehensive discussion of 
renormalons see [M. Beneke, Phys. Rep. 317, f (f999)]), and (ii) a contour-improved combination 
of the Borel method with an analytic continuation by conformal mapping, and Pade approximations 
in the conformal variable. The singularity structure in the case of the LoSurdo-Stark effect in 
the complex Borel plane is shown to be similar to (divergent) perturbative expansions in quantum 
chromodynamics. 
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I. INTRODUCTION 

Consider the energy shift of the ground state of atomic 
hydrogen in an electric field of field strength F which 
we assume to lie along the z-axis: the energy displace- 
ment can be expressed in perturbation theory as a for- 
mal power series. The first nonvanishing perturbation [in 
atomic units, see also Eq. (|9|) below] is the second-order 
effect 
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where the sum over m runs over the entire spectrum, 
including the continuum, but excluding the IS ground 
state, and E m is the nonrelativistic (Schrodinger) energy 
of the mth state. A well known, but perhaps surprising 
result says that the coefficients of the terms of order F 4 , 
F 6 , F 8 , . . . grow so rapidly that the series in F ultimately 
diverges, irrespective of how small the field strength is. 
The convergence radius of the factorially divergent per- 
turbation series is zero. The resummation of the diver- 
gent series is problematic in the considered case, because 
the Borel transform, from which the physically correct, 
finite result is obtained by evaluating the Laplace-Borel 
integral - see Eq. ( [ll|) in Sec. Ill below -, exhibits a rich 
singularity structure in the complex plane. 

The purpose of this paper is to present numerical ev- 
idence that divergent perturbation series whose Borel 
transform exhibits a rich singularity structure in the com- 
plex plane, can be resummed to the complete, physi- 
cally relevant result. The resummation methods use as 
input data only a finite number of perturbative coeffi- 
cients. Problematic singularities on the positive real axis 



in the Borel plane are treated by appropriate integra- 
tion prescriptions, fn particular, it is shown that the 
Borel transform of the divergent perturbation series for 
the LoSurdo-Stark effect involves two cuts in the Borel 
plane, generated essentially by the divergent alternating 
and nonalternating subcomponents of the perturbation 
series. This singularity structure is also expected of the 
(divergent) perturbation series in quantum field theory, 
notably quantum chromodynamics (in this case, the al- 
ternating and nonalternating factorially divergent com- 
ponents correspond in their mathematical structure to 
the ultraviolet and infrared renormalons). 

We present results which suggest that the integration 
contours and resummation techniques discussed here may 
be of relevance, at least in part, to theories with degener- 
ate minima. As a byproduct of these investigations, nu- 
merical pseudo-eigenvalues are obtained for the LoSurdo- 
Stark effect; selected field strengths and atomic states are 
considered. 

The LoSurdo-Stark effect and its associated divergent 
perturbative expansion including the the nonperturba- 
tive, nonanalytic imaginary contributions have attracted 
considerable attention, both theoretically and experi- 
mentally fl|-|34||. Experiments have been performed in 
field strengths up to a couple of MV/cm p5|-p8[. One 
might be tempted to say that the autoionization decay 
width could be interpreted as a paradigmatic example for 
a nonperturbative effect which exhibits fundamental lim- 
itations to the validity of perturbation theory (unless the 
perturbative expansion is combined with resummation 
methods). We briefly summarize here: The Rayleigh- 
Schrodinger perturbation series for the LoSurdo-Stark 
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effect |39|, [40[ can be formulated to arbitrarily high or- 
der Q]. The perturbative coefficients grow factorially in 
absolute magnitude M, and the radius of convergence 
of the perturbation series about the origin is zero. The 
perturbation series is a divergent, asymptotic expansion 
in the electric field strength F, i.e. about zero electric 
field. This means that the perturbative terms at small 
coupling first decrease in absolute magnitude up to some 
minimal term. After passing through the minimal term, 
the perturbative terms increase again in magnitude, and 
the series ultimately diverges. 

By the use of a resummation method, it is possible to 
assign a finite value to an otherwise divergent series, and 
various applications of resummation methods in math- 
ematics and physics have been given, e.g., in jll] [tHj . 
When a divergent series is resummed, the superficial pre- 
cision limit set by the minimal term can be overcome, 
and more accurate results can be obtained as compared 
to the optimal truncation of the perturbation series (see 
also the numerical results in the Tables of The di- 

vergent perturbation series of the LoSurdo-Stark effect 
has both alternating and nonalternating components (as 
explained in Sec. O below). The resummation of non- 
alternating series or of series which have a leading or 
subleading divergent nonalternating component, corre- 
sponds to a resummation "on the cut" in the complex 
plane (4l],|42|. 

Rather mathematically motivated investigations re- 
garding the Borel summability of the divergent perturba- 
tion series for the LoSurdo-Stark effect were performed 
in Q and jlTj , and it was established that the perturba- 
tion series of the LoSurdo-Stark effect is Borel summable 
in the distributional sense (for the definition of "distri- 
butional Borel summability" we refer to [Q). Here, to 
supplement the mathematically motivated investigations, 
we consider the calculation of transforms of the diver- 
gent series, which use as input data only a finite number 
of perturbative coefficients and exhibit apparent conver- 
gence to the complete, physically relevant result. 

In the remarkable investigation |^5| , whose significance 
may not have been sufficiently noticed in the field of 
large-order perturbation theory, it was not only shown 
that it is possible to perform the required analytic con- 
tinuation of the Borel transform beyond its radius of con- 
vergence by employing Pade approximants, but that it 
is also possible to reconstruct the full physical result, 
including the imaginary contribution which corresponds 
to the autoionization decay width, by integration in the 
complex plane. 

In Sec. |l[ we discuss the singularity structure of the 
Borel transform in the complex plane. The structure of a 
doubly-cut plane has been postulated for quantum chro- 
modynamic perturbation series fl49|, |5(i|] , and this struc- 
ture has been exploited to devise resummation prescrip- 
tions based on conformal mappings p9|-p5[. Here, we 
present results which suggest that the convergence of the 
transforms obtained by conformal mapping can be im- 
proved if Pade approximants in the conformal variable 



are used (see also |5f|). We also discuss improvements 
of the "pure" Borel-Pade method (these additional im- 
provements take advantage of leading renormalon poles). 
Also, in comparison to the investigation (29], we use here 
a slightly modified, but equivalent integration contour for 
the evaluation of the generalized Borel integral (see fiq| 
and Sec. [HI] below). Our version of the integration con- 
tour exhibits the additional terms which have to be added 
to the otherwise recommended principal-value prescrip- 

tion mum. 

As stressed above, it has been another main motivation 
for the current investigation to establish the singularity 
structure of the Borel transform in the complex plane, 
and to demonstrate the analogy of the singularity struc- 
ture of the perturbation series for the LoSurdo-Stark ef- 
fect to quantum chromodynamic perturbation series. We 
also consider a divergent perturbation series generated by 
a model problem for theories with degenerate minima. In 
the particular model case discussed here, a perturbation 
series with real coefficients is summed to a real result 
- in contrast to the LoSurdo-Stark effect, there is no 
imaginary part involved in this case. One of the three 
alternative integration contours introduced in has to 
be employed. 

This paper is organized as follows: In Sec. E3, we 
give a brief outline of the perturbative expansion for the 
LoSurdo-Stark effect. In Sees. |n] and |v|, we describe 
the resummation methods which are used to obtain the 
numerical results presented in Sec. In Sec. |vi| , we 
consider theories with degenerate mini ma. We conclude 
with a summary of the results in Sec. |VIl| . Finally, the 
connection of the current investigation to quantum field 
theoretic perturbation series and to double-well oscilla- 
tors are discussed in the Appendixes ^ and [B|. 



II. PERTURBATION SERIES 
FOR THE LoSURDO-STARK EFFECT 

In the presence of an electric field, the SO(4) symmetry 
of the hydrogen atom is broken, and parabolic quantum 
numbers n\, and m are used for the classification of 
the atomic states |57j]. For the Stark effect, the pertur- 
bative expansion of the energy eigenvalue E(n\ , ni , m, F) 
reads [see Eq. (59) of §], 



E( ni ,n 2 ,m,F) ~ Ei% m F N , 



(1) 



where F is the electric field strength. For N — > oo, the 
leading large-order factorial asymptotics of the perturba- 
tive coefficients have been derived in |L6| as 

£?W +(-i)N A (N) N^oo (2) 



where A. 



(N) 
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is given as an asymptotic series, 
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x J2al injTn (2n j + m + N-k)\. (3) 

fc=0 

The quantities a ^™ 3?rl are constants. The if-coefficients 
in Eq. (^|) are given by 

K(rii, rij, rn, N) = — [27m 3 n^! (nj + to)!] 

xe 3 cp{3(n i -n i )} 6 2 ^ +m+1 (3n 3 /2) JV . (4) 

Here, the principal quantum number n as a function of 
the parabolic quantum numbers n\, n 2 and m is given 
by [see Eq. (65) in §] 
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FIG. 1. Integration contour C+i for the evaluation of 
the generalized Borel integral defined in Eq. Poles 
displaced from the real axis are evaluated as full poles, 
whereas those poles which lie on the real axis are treated 
as half poles. 



ni + n 2 + |to| + 1 . 



(5) 



According to Eq. (0), the perturbative coefficients 



En^l 2 m, for large order N — > oo of perturbation theory, 
can be written as a sum of a nonalternating factorially di- 
vergent series [first term in Eq. (Q)] and of an alternating 
factorially divergent series [second term in Eq. (Q)]. Be- 
cause the a™*™ 3 " 1 i n Eq, (fj) arc multiplied by the factorial 
(2rij + to + N — A;)!, we infer that for large perturbation 
theory order N, the term related to the aQ ,rljm coefficient 
(k = 0) dominates. Terms with k > 1 are suppressed in 
relation to the leading term by a relative factor of 1 /N k 
according to the asymptotics 



(2nj + m + N - k)\ 1 
A* 7 
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■ N)\ 



(0) 



for A — *■ oo. The leading (fc = 0) -coefficient has been 
evaluated in M as 



O n = 1 



(7) 



According to Eqs. (||), (||) and (0), for states with m < 
n 2 , the nonalternating component of the perturbation 
series dominates in large order of perturbation theory, 
whereas for states with ri\ > n 2l the alternating compo- 
nent is dominant as A — ► oo. For states with m = n 2 , 
the odd- TV perturbative coefficients vanish fl6[| , and the 
even- A coefficients necessarily have the same sign in large 
order [see Eq. (§)]. According to Eq. (g), there are sub- 
leading divergent nonalternating contributions for states 
with ri\ > ri2, and there exist subleading divergent al- 
ternating contributions for states with n\ < n 2 . This 
complicates the resummation of the perturbation series. 



III. BOREL-PADE RESUMMATION 

The resummation of the perturbation series ([!]) by a 
combination of the Borel transformation and Pade ap- 
proximants proceeds as follows. First we define the pa- 
rameter 



A = 2 max(ni , n 2 ) 



1 



(8) 



The large-order growth of the perturbative coefficients 
[see Eqs. (||) and (||)] suggests the definition of the (gen- 
eralized) Borel transform [see Eq. (4) in |5^] j 

E B (z) = E B (ni , n 2 , to, z) 



= B {1X) [J5(m,n a ,m); z] 



E 

N=0 



(N) 

nmztn TV 



r(A + A) 



(9) 



where we consider the argument z of E B (z) as a com- 
plex variable and A is defined in Eq. (||). The additive 
constant (in this case A) in the argument of the Gamma 
function is chosen in accordance with the notion of an 
"asymptotically improved" resummation (see also p8|). 
ft is observed that the additive constant A can be shifted 
by a small integer without affecting the convergence of 
the Borel resummed series. Because the perturbative co- 



sfficients E. 



(N) 



diverge factorially in absolute magni- 



tude, the Borel transform Eq (z) has a finite radius of con- 
vergence about the origin. The evaluation of the (gener- 
alized) Laplace-Borel integral [see Eq. ( |ll| ) below] there- 
fore requires an analytic continuation of E B (z) beyond 
the radius of convergence. The "classical" Borel integral 
is performed in the z-range z € (0, oo), i.e. along the pos- 
itive real axis [see e.g. Eqs. (8.2.3) and (8.2.4) of Q]. ft 
has been suggested in f5q] that the analytic continuation 
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of ([|) into regions where F retains a nonvanishing, albeit 
infinitesimal, imaginary part can be achieved by evalu- 
ating Pade approximants. Using the first M + 1 terms 
in the power expansion of the Borel transform Eb(z), we 
construct the Pade approximant (we follow the notation 
of gj) 



V M (z) 



[M/2] / \{M + l)/2] 



(*)> 



(10) 



E B 



where [x] denotes the largest positive integer smaller 
than x. We then evaluate the (modified) Borel integral 
along the integration contour C+i shown in Fig. |l| in or- 
der to construct the transform TEm (F) where 



TEm(F) 



dtt^ 1 exp(~t)V M {Ft) . 



(11) 



The successive evaluation of transforms TEm{F) in in- 
creasing transformation order M is performed, and the 
apparent convergence of the transforms is examined. 
This procedure is illustrated in Tables 1 and II of p6[ . 
In the current evaluation, a slightly modified scheme is 
used for selecting the poles in the upper right quadrant 
of the complex plane as compared to Q . 

The contour C+i is supposed to encircle all poles at 
t = Zi in the upper right quadrant of the complex plane 
with arg zi < 7r/4 in the mathematically negative sense. 
That is to say, the contribution of all poles z% with 
He Zi > 0, Imzi > and Rez^ > Imzj, 



-2 ?ri 



Res t 



A-l 



exp(-t) V M (Ft) 



is added to the principal value (P.V.) of the integral (|l] 
carried out in the range t € (0, oo). Note the further 
restriction (Imz^ < Rez^ or equivalently arg z% < tt/4) 
with regard to the selection of poles in comparison to 
the previous investigation p6| . In practical calculations, 
this modification is observed not to affect the numerical 
values of the transforms TEm(F) defined in Eq. ([ll]) in 
higher transformation order M > 10 [i.e. for large M, see 
also Eq. ( p^ ) below], because the poles are observed to 
cluster near the real axis in higher transformation order, 
and so the contribution of poles with 7r/4 < argz; < it/2 
gradually vanishes. We have, 

/•OO 

TE M {F) = (P.V.) / dtt^ 1 exp(~t)V M {Ft) 
Jo 

-2 7ri V Resi A_1 exp(-t)V M (Ft) . (12) 



The principal- value prescription [first term in Eq. (12)] 
for the evaluation of the Laplace-Borel integral has been 
recommended in |5^,|5^]. This prescription leads to a 
real (rather than complex) result for the energy shift 
and cannot account for the width of the quasistationary 
state. The additional pole contributions [second term in 



Eq. ([12])] are responsible for complex-valued (imaginary) 
corrections which lead, in particular, to the decay width. 

By contour integration (Cauchy Theorem) and Jor- 
dan's Lemma, one can show that the result obtained 
along C+i is equivalent to an integration along the 
straight line with argz = tt/4, 



/>oo 

TE M (F) = c x / dtt*- 1 exp(-ct)V M (F 
Jo 



ct) (13) 



where c = exp(i 7r/4). This contour has been used in |2q] 
(see also p. 815 in p0f). The exponential factor exp(— ct) 
and the asymptotic behavior of the Pade approximant 
Vm{F ' ct) as t — ► 00 together ensure that the integrand 
falls off sufficiently rapidly so that the Cauchy Theorem 
and Jordan's Lemma can be applied to show the equiva- 
lence of the representations ( |l2] ) and ([l3]) . 

The representation (|ll]) illustrates the fact that the 
integration in the complex plane along C+i analytically 
continues the resummed result in those cases where the 
evaluation of the standard Laplace-Borel integral is not 
feasible due to poles on the real axis. The representations 
( pi] ) and ( p"2| ) serve to clarify the role of the additional 
terms which have to be added to the result obtained by 
the principal-value prescription in order to obtain the 
full physical result, including the nonperturbative, non- 
analytic contributions. Note that, as stressed in fl46| , 
the pole contributions in general do not only modify the 
imaginary, but also the real part of the resummed value 
for the perturbation series. 

Formally, the limit of the sequence of the TEm{F) as 
M — * 00, provided it exists, yields the nonperturbative 
result inferred from the perturbative expansion (hi), 



M 



lim TEm (F) = E(F) = E(ni,n 2 ,m,F) 



(14) 



Because the contour C+i shown in Fig. [j] extends into 
the complex plane, the transforms TEm(F) acquire an 
imaginary part even though the perturbative coefficients 
in Eq. (Q) are real. 

In the context of numerical analysis, the concept of 
incredulity |3l| may be used for the analysis of the con- 
vergence of the transforms TEm(F) of increasing order 
M. If a certain number of subsequent transforms ex- 
hibit apparent numerical convergence within a specified 
relative accuracy, then the calculation of transforms is 
stopped, and the result of the last calculated transfor- 
mation is taken as the numerical limit of the series under 
investigation. It has been observed in [[l6],[5(| that for 
a number of physically relevant perturbation series, the 
apparent numerical convergence of the transforms ([ll]), 
with increasing transformation order, leads to the phys- 
ically correct results. 

It is observed that the rate of convergence of the trans- 
forms (fill) can be enhanced if instead of the unmodified 



Pade approximants (10) leading "renormalon" poles are 
explicit used for the construction of modified approxi- 
mants. For the ground state, this entails the following 
replacement in Eq. (lU 
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V M {z) -» V'm(z) 



where 
V'm(z) 



l-z 2 



r M + 4 T 



r M-3, 



(*), (15) 



where -Eb(C) = (1 — C 2 ) ^b(C)- F° r the excited state with 
quantum numbers ri\ = 3, rii = and m = 1, we replace 



V M {z) -» P"m(z) . 



where 



1 



(z) , (16) 



E" 



where Eq(Q = (1 — ( 2 ) Eb(0- The resummation method 
by a combination of Borel and Pade techniques - current 
Section - will be referred to as "method I" throughout 
the current paper. 



IV. DOUBLY-CUT BOREL PLANE 
AND RESUMMATION 
BY CONFORMAL MAPPING 

According to Eqs. (0) and (0), the perturbative coeffi- 



cient E, 



(JV) 



for large N, can be written as the sum of 
an alternating and of a nonalternating divergent series. 
In view of Eqs. M) and (0), we conclude that the series 
defined in Eq. (g)^ 



Eb(z) = 



E. 



(N) 



n i n : > at. _/V 



N=0 

has a radius of convergence 



^ T(N + A) 



2 
3^3 



(17) 



about the origin, where n is the principal quantum num- 
ber [see Eq. (H)]. Therefore, the function 



£ B (w) 



E y 



(N) 



^ r(iv + A) 



w 



N 



(18) 



has a unit radius of convergence about the origin. It is not 
a priori obvious if the points w = — 1 and w = +1 repre- 
sent isolated singularities or branch points. The asymp- 
totic properties (Q) and (||) together with Eq. (||) suggest 
that the points w = — 1 and w = +1 do not constitute 
poles of finite order. We observe that the leading facto- 
rial growth of the perturbative coefficients in large per- 
turbation order N is divided out in the Borel transform 
(|l8|), which is a sum over N. The perturbative coeffi- 
cient E^n 2 m can be written as an asymptotic series over 



k [see Eq. (0)]. We interchange the order of the summa- 
tions over N and k, we use Eq. (||) and take advantage 
of the identity 



oo N 

E w = Li *H- 



JV=0 



(19) 



The Borel transform £b(w) can then be written as a sum 
over terms of the form Tk(w) where for k — > oo, 



T fc (u;) ~ C(ni,rij,m) a k % 3 Li k (w) . 
The coefficient C (rij , nj , to) is given by 

C (rij , rij , to) = — [27rn 3 nj ! (to,- + to) !] 
xexp {Zim-nj)} 6 2 ^+ m+1 . 



(20) 



(21) 
-1 



These considerations suggest that the points w - 
and w = +1 represent essential singularities (in this case, 
branch points) of the Borel transform £b{w) defined in 
Eq. (p"8|). For the analytic continuation of £b{w) by con- 
formal mapping, we write w as 



2y 



1 



(22) 



(this conformal mapping preserves the origin of the com- 
plex plane). Here, we refer to w as the Borel variable, 
and we call y the conformal variable. We then express 
the Mth partial sum of (pit) as 



M 



= E 



(TV) 



N=0 



T(N + A) 



w 



N 



M 

E 

N=0 



C N y N +0(y M+ ') 



(23) 



where the coefficients Cm are uniquely determined [see, 
e.g., Eqs. (36) and (37) of @]. We define the partial 
sum of the Borel transform, re-expanded in terms of the 
conformal variable y, as 



M 



£ B M (y) = J2 c *y N 



(24) 



N=0 



We then evaluate (lower-diagonal) Pade approximants to 
the function ^(y), 



cllM 
C B 



(y) 



[M/2]/p/ + l)/2] 



(?/) 



(25) 



We define the following transforms, 

T"E M (F) = s x f dw w^ 1 exp(-w) £b M (y(w)) 
Jc +1 



(26) 
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At increasing M, the limit as M — > oo, provided it ex- 
ists, is then again assumed to represent the complete, 
physically relevant solution, 



E(F) = lim T"E M (F). 

M — > oo 



(27) 



We do not consider the question of the existence of this 
limit here (for an outline of questions related to these 
issues we refer to |50[ ; potential problems at excessively 
strong coupling are discussed in Sec. IIC of J62[ ). 
Inverting Eq. © yields [see Eq. ©] 



y{w) 



-VI 



VT+^+\/r 



(28) 



The conformal mapping given by Eqs. (|22|) and ( |28| ) maps 
the doubly cut w-plane with cuts running from w = 1 to 

= oo and w = -1 to w = -oo unto the unit circle in 
the complex y-plane. The cuts themselves are mapped 
to the edge of the unit circle in the y-plane. 

In comparison to the investigations Q and fl50| , 
we use here a different conformal mapping defined in 
Eqs. (p||) and (|2^) which reflects the different singular- 
ity structure in the complex plane [cf. Eq. (27) in p9|]. 
We also mention the application of Pade approximants 
for the numerical improvement of the conformal map- 
ping performed according to Eq. (pq). In comparison to 
a recent investigation |5q| , where the additional Padc- 
improvement in the conformal variable is also used, we 
perform here the analytic continuation by a mapping 
whose structure reflects the double cuts suggested by 
the asymptotic properties of the perturbative coefficients 
given in Eqs. @, (|) and (§) [cf. Eq. (5) in 

The method introduced in this Section will be referred 
to as "method II". It is one of the motivations for the 
current investigation to contrast and compare the two 
methods I and II. A comparison of different approaches 
to the resummation problem for series with both alter- 
nating and nonalternating divergent components appears 
useful, in part because the conformal mapping (without 
further Pade improvement) has been recommended for 
the resummation of quantum chromodynamic perturba- 
tion scries {|||5(|. 

We do not consider order-dependent mappings 
here [|5l|-[5"4| . For an order-dependent mapping to be con- 
structed, the conformal mapping in Eq. ( |22] ) has to be 
modified, and a free parameter p has to be introduced. 
The coefficients Cn in the accordingly modified Eq. (|24| ) 
then become p-dependent. The free parameter p is cho- 
sen so that the p-dependent coefficient Cm (p) of order M 
vanishes. Consequently, the choice of p depends on the 
order M of perturbation theory, and in this way the map- 
ping becomes order-dependent. Certain complications 
arise because p cannot be chosen arbitrarily, but has to 
be selected, roughly speaking, as the first zero of the p- 
dependent coefficient Cm (p) for which the absolute mag- 
nitude of the derivative C' M (p) is small (this is explained 
in M, p. 886). It is conceivable that with a judicious 



choice of p, further acceleration of the convergence can 
be achieved, especially when the order-dependent map- 
ping is combined with a Pade approximation as it is done 
here in Eq. (^5|) with our order-independent mapping. In 
the current investigation, we restrict the discussion to 
the conformal order-independent mapping ( |22j ) which is 
nevertheless optimal in the sense discussed in pS.OT. 



V. NUMERICAL CALCULATIONS 

In this section, the numerical results based on the re- 
summation methods introduced in Sees. Ill and IV are 
presented. 



Before we describe the calculation in detail, 
we should note that relativistic corrections to both the 
real and the imaginary part of the energy contribute at 
a relative order of (Za) 2 compared to the leading non- 
relativistic effect which is treated in the current inves- 
tigation (and in the previous work on the subject, see 
e.g. |l(],gf|). Therefore, the theoretical uncertainty due 
to relativistic effects can be estimated to be, at best, 1 
part in 10 4 (for an outline of the relativistic and quan- 
tum electrodynamic corrections in hydrogen see [|63|-|69| ) . 
Measurements in very high fields are difficult |35| . At the 
achievable field strengths to date (less than 0.001 a.u. or 
about 5MV/cm), the accuracy of the theoretical predic- 
tion exceeds the experimental precision, and relativistic 
effects do not need to be taken into account. 

The perturbative coefficients En^l 2 m defined in Eq. (^) 
for the energy shift can be inferred, to arbitrarily high or- 
der, from the Eqs. (9), (13-15), (28-33), (59-67) and (73) 
in Q. The atomic unit system is used in the sequel, as 
is customary for this type of calculation [|], ||, ^|, |ll| . The 



unit of energy is a 



27.211 eV where a is the 



fine structure constant, and the unit of the electric field 
is the field strength felt by an electron at a distance of 
one Bohr radius aeohr to a nucleus of elementary charge, 
which is l/(47re ) (e/a| olir ) = 5.142 x 10 3 MV/cm (here, 
eo is the permittivity of the vacuum). 

We consider the resummation of the divergent pertur- 
bative expansion ([!]) for two states of atomic hydrogen. 
These are the ground state (n\ = = m = 0) and an 
excited state with parabolic quantum numbers ni = 3, 
U2 = 0, m = 1. We list here the first few perturbative 
coefficients for the states under investigation. For the 
ground state, we have (in atomic units), 



Eqoo{F) - -- - - F - 



1 9 „ 2 3 555 4 



2 512 779 
512 



F 



64 

13 012 777 803 
16 384 



F 



(29) 



The perturbation series for the state n\ = 3, ri2 = 0, 
m = 1 is alternating, but has a subleading nonalternat- 
ing component [see Eq. (^J)]. The first perturbative terms 
read 



6 



31875 



F 



54 140 625 



715 751953125 



1G 



F 



(30) 



Note that for F = 0, the unperturbed nonrelativistic en- 
ergy is recovered, which is — 1/(2 n 2 ) in atomic units. In 
contrast to the real perturbative coefficients, the energy 
pseudoeigenvalue (resonance) E(m, ri2, m, F) has a real 
and an imaginary component, 

2, 

E(nx,n 2 ,m,F) = ReE nin2Tn (F) - - T niTl2m (F) , (31) 

where T nin2m (F) is the autoionization width. 

Using a computer algebra system |70| , |7l| , the first 50 
nonvanishing perturbative coefficients are evaluated for 
the ground state, and for the state with parabolic quan- 
tum numbers rii = 3, n 2 = 0, m = 1, we evaluate the 
first 70 nonvanishing perturbative coefficients. The ap- 
parent convergence of the transforms defined in Eqs. ( pd| ) 
and (53) in higher order is examined. In the case of the 
Borel-Pade transforms defined in Eq. (pi]), use is made 
of the replacements in Eqs. ( |l5| ) and (|l^) pleading renor- 
malon poles are being put in by hand"]. This procedure 
leads to the numerical results listed in Tables [j] and ||. 
The numerical error of our results is estimated on the 
basis of the highest and lowest value of the four highest- 
order transforms. 

An important result of the comparison of the meth- 
ods introduced in Sees. Ill and IV is the following: Both 
methods appear to accomplish a resummation of the per- 
turbation series to the physically correct result. Method I 
(Borel+Pade with leading renormalon poles, see Sec. Ill) 
and method II (Borel+Pade-improved conformal map- 



ping, see Sec. IV) appear to lead to results of comparable 
accuracy. 

To date, a rigorous theory of the performance of the 
resummation methods for divergent series of the type dis- 
cussed in this work (with alternating and nonalternating 
components) does not exist. The logarithmic singularities 
introduced by the branch points of higher-order polylog- 
arithms [see the index k in Eq. (19)] are difficult to ap- 
proximate with the rational functions employed in the 
construction of Pade approximants. A solution to the 
problem of approximating the logarithmic singularities, 
based on finite number of perturbative coefficients, would 
probably lead to further optimizimation of the rate of 
convergence of the transformed series. Within the cur- 
rent scheme of evaluation, the problematic logarithmic 
singularities may be responsable, at least in part, for cer- 
tain numerical instabilities at higher transformation or- 
der, e.g. in the result for T"E 70 {F = 2.1393 x 10~ 4 ) in 
Eq. (|32|) below. 

For the atomic state with quantum numbers n\ = 3, 
rii = and m = 1, the evaluation of the transforms 
TEm{F) defined in Eq. ( |TT| ) (method I) and of the trans- 
forms T"Em(F) defined in Eq. (^6|) (method II) in trans- 
formation order M — 67, 68, 69, 70 for a field strength of 



F = 2.1393 x 10~ 4 . Method I leads to the following 
results, 

TE 67 (F = 2.1393 x 10~ 4 ) = -0.015 860 468199 2 
-i 0.529 048 x 10~ 6 , 

TE 6S (F = 2.1393 x 10~ 4 ) = -0.015 860 468 200 9 
-i 0.529 047 x 10~ 6 , 

TE 6g (F = 2.1393 x 10~ 4 ) = -0.015 860 468 198 9 
-i 0.529 048 x 10~ 6 and 



TE 70 (F = 2.1393 x 10~ 4 ) = -0.015 860 468194 5 



-i 0.529 015 x 10 



-r> 



(32) 



TABLE I. Real and imaginary part of the energy pseu- 
doeigenvalue Eooo(F) for the ground state of atomic hydrogen 
(parabolic quantum numbers rii = 0, ni = 0, m = 0). 



F (a.u.) 


ReSoooCF) 


roooCF) 


0.04 


-0.503 771 591 013 654 2(5) 


3.892 699 990(1) x 10 


0.06 


-0.509 203 451 088(2) 


5.150 775 0(5) x 10~ 4 


0.08 


-0.517 560 50(5) 


4.539 63(5) X 10~ 3 


0.10 


-0.527 419 3(5) 


1.453 8(5) x 10~ 2 


0.12 


-0.537 334(5) 


2.992 7(5) x 10~ 2 


0.16 


-0.555 24(5) 


7.131(5) x 10~ 2 


0.20 


-0.570 3(5) 


1.212(5) x 10 _1 


0.24 


-0.582 6(1) 


1.767(5) x 10 _1 


0.28 


-0.591 7(5) 


2.32(3) x 10" 1 


0.32 


-0.600(5) 


2.92(3) x 10 _1 


0.36 


-0.604(5) 


3.46(3) x 10 _1 


0.40 


-0.608(5) 


4.00(5) x 10 _1 
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Method II yields the following data, 
T"E 67 (F = 2.1393 x 1CT 4 ) = -0.015 860 468 200 4 
-i 0.529 047 x 10~ 6 , 

T"E 6S {F = 2.1393 x 10" 4 ) = -0.015 860 468 200 3 
-i 0.529 047 x 10~ 6 , 



T"E 69 (F = 2.1393 x 10" 



-0.015 860 468 200 4 



-i 0.529 047 x 10~ 6 and 

T"E 70 {F = 2.1393 x 10" 4 ) = -0.015 860 468 203 3 

-i 0.529 046 x 10~ 6 . (33) 

Numerical results obtained by resummation are pre- 
sented in Tables | and || for a variety of field strengths 
and for the two atomic states under investigation here. 
Results are compared to the numerical calculation |5). 
which yields very accurate data for all experimentally 
accessible electric field strengths to date. In addition, 
it should be noticed that the inaccuracies at excessively 
large field of have been pointed out by the same au- 
thors in ^7j. However, not all atomic states considered 
in 11 were treated in the later investigation ^7|. Our 



are con- 
How- 
7j leaves 
3, 



27 1 



data for the ground state indicated in Table 
sistent with the numerical results obtained in 
ever, it should be noted that the later work 
out the excited state with quantum numbers n\ 
112 = and m — 1 for which results are given here here 
in Table |l[ To the best of our knowledge, the numeri- 
cal discrepancy with || for the excited state with quan- 
tum numbers ri\ = 3, rii = and m = 1 has not been 
recorded in the literature. We do not claim here that 
it would have been impossible to discern this discrep- 
ancy with the other methods which have been devised 
for the theoretical LoSurdo-Stark problem. Notably, it 
appears likely that the approach from |27| or the method 
presented in could easily be generalized to the par- 
ticular excited state considered here, and that such a 
generalization would lead to very accurate results. We 
merely include Table |l| here in order to illustrate the util- 
ity of the rather unconventional resummation method for 
the regime of large coupling, where even rather sophisti- 
cated numerical methods, which avoid the intricacies of 
a small-field perturbative expansion, have been shown to 
be problematic_0, . We confirm that the numerical 
data given in || are accurate up to a field strength of 
F sa 0.1 for the ground state and up to F w 3 x 10~ 4 
for the excited (n = 5)-state with rii = 3, n-z = and 
m = 1. 



TABLE II. Real part and imaginary part of the energy pseudo-eigenvalue Ezo\{F) for the excited state with parabolic 
quantum numbers n\ = 3, ri2 = 0, m = 1. The field strength F is given in atomic units. The data are compared to Ref. 
Discrepancies are observed at large electric field strength. 









Real part of the resonance Rei?3oi(-F) 


Autoionization 


decay width r3oi(-F) 


F (a.u.) 






Ref. |] 


Our results 


Ref. | 


Our results 


1.5560 x 


10" 


-4 


-0.016 855 237 2 


-0.016 855 237 140 761 7(5) 


0.42 x 10~ 9 


0.421 683(5) x 10 -9 


1.9448 X 


10" 


-4 


-0.016 179 388 5 


-0.016 179 388 257 0(5) 


0.143 8 x 10 -6 


0.143 773(5) x 10 -6 


2.1393 x 


10" 


-4 


-0.015 860 468 


-0.015 860 468 20(1) 


0.105 7 x 10~ 5 


0.105 09(5) x 10~ 5 


2.5282 x 


10" 


-4 


-0.015 269 204 


-0.015 269 293(1) 


0.175 60 x 10~ 4 


0.176 39(5) x 10~ 4 


2.9172 x 


10" 


-4 


-0.014 740 243 


-0.014 742 60(3) 


0.976 51 x 10~ 4 


0.999 96(9) x 10~ 4 


3.3061 x 


10" 


-4 


-0.014 242 49 


-0.014 260 2(3) 


0.278 53 x 10~ 3 


0.295 4(2) x 10~ 3 
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TABLE III. Resummation of the asymptotic series 
for the generating functional of a zero-dimensional the- 
ory with degenerate minima given in Eqs. ( |38| ) and (^9|). 
We have g — 0.01. Results in the third column are ob- 
tained by the method indicated in Eq. (Q) along the 
integration contour Co (see (46)). The partial sums in 
the second column are obtained from the asymptotic se- 
ries @. 



M 


partial sum 


TZ M ( 


9 = 0.01) 


2 


1.081 000 


1.102 


326 


3 


1.094 860 


1.096 


141 


4 


1.108 373 


1.089 


875 


5 


1.125 832 


1.090 


695 


6 


1.153 942 


1.092 


000 


7 


1.208 154 


1.091 


596 


G 
O 




1.091 


389 


9 


1.642 718 


1.091 


553 


10 


2.545 239 


1.091 


545 


11 


5.438 230 


1.091 


503 


12 


1.5 x 10 1 


1.091 


525 


13 


5.5 x 10 1 


1.091 


527 


14 


2.2 x 10 2 


1.091 


519 


15 


9.5 x 10 2 


1.091 


523 


16 


4.5 x 10 3 


1.091 


523 


17 


2.2 x 10 4 


1.091 


521 


18 


1.2 x 10 5 


1.091 


522 


19 


6.9 x 10 5 


1.091 


522 


20 


4.1 x 10 6 


1.091 


522 


exact 


1.091 522 


1.091 522 



VI. MODEL EXAMPLE FOR DEGENERATE 
MINIMA 



r(2iV + 1/2) 

r(jv + 1) 



T(N) 



1 + 0{ N 



(36) 



it is easy to explicitly establish the factorial divergence of 
the series (see also p. 888 of Q). The generating func- 
tional in zero dimensions has been proposed as a paradig- 
matic example for the divergence of perturbation theory 
in higher order. It can be resummed easily to the non- 
perturbative result; in particular it is manifestly Borcl 
summable, and no singularities are present on the posi- 
tive real axis. 

Complications are introduced by degenerate minima. 
As a second example, we consider the modified generat- 
ing functional [compare with Eq. (2.6) on p. 15 of |73j | 
and with Eq. (40.1) on p. 854 of §§]: 



OO 

f d$ 



exp 



(37) 



The expansion of the exponential in powers of the cou- 
pling g leads to a divergent asymptotic series, 



{ ' AT! J V2^ V y 2 7 



E 

N=0 



N 



d$ 



-1/2* 2 



N 

Er(2iv-j + i) 



2N-j\ t ^ , 3 x2(jv-i) fgd*^ ' 



We consider the generating functional in a zero- 
dimensional theory (in this case, the usual path integral 
reduces to an ordinary integral). First, we briefly con- 
sider the $ 4 -theory in zero dimensions [see Eq. (9-177) 
ff. in 173]]; the generating functional reads 



d$ 

V5n 



exp 



5$ 



(34) 



The strictly alternating divergent asymptotic expansion 
in powers of g for g — > reads, 



A N T{2N + 1/2) 



E 



(-.9) 



N 



(35) 



On using the known asymptotics valid for N — ► oo, which 
is this case yield the "large-order" asymptotics of the per- 
turbative coefficients, 



J2 2^ { J } g N 



N=0 



T(N - 1/2) 



oo N 



E 

N=0 



8 N T(2N+l/2) N 
/nT(N + l) 9 



(38) 



where C^ [ {x) denotes a Gegenbauer (ultraspherical) 
polynomial. Note that terms of half-integer power of g 
entail an odd power of the field and vanish after integra- 
tion. The first few terms of the asymptotic expansion 
read, 

Z'($) = 1 + 6 g + 210 g 2 + 13860 g 3 
+ 1351350 g 4 + 174594420 g 5 
+28109701620 g 6 + 5421156741000 g 7 
+ 1218404977539750 g s + 0{g 9 ). (39) 
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For the perturbative coefficients 



Cv — 



8" T(2N+ 1/2) 



we establish the following asymptotics, 



C 



1 



N 



N-'-32 N T(N + 1) . 



ttV2 



(40) 



(41) 



Due to the nonalternating character of the expansion 
(|38|), it is not Borel summable in the ordinary sense. 
Rather, it is Borel summable in the distributional 
sense [[ItJQ. Here, we present numerical evidence sup- 
porting the summability of the divergent expansion (|39| ) 
based on a finite number of perturbative coefficients. The 
final integration is carried out along the contour Co in- 
troduced in |4(J [see also Eq. ( p4| ) below]. The same 
contour has also been used for the resummation of diver- 
gent perturbation series describing renormalization group 
(anomalous dimension) 7 functions ]5q] . As explained 
in ^6|, the integration along Co, which is based on the 
mean value of the results obtained above and below the 
real axis, leads to a real final result if all perturbative 
coefficients are real. 

In particular, the resummation of the divergent ex- 
pansion ( ^ ) is accomplished as follows. We first define 
the Borel transform of the generating functional by [see 
Eq. (4) in Q and the discussion after Eq. M)] 



Z' B {Z) EE SC 1 ' 1 ) [Z >. z] 



E 



^ o r(iv + i) 



(42) 



Pade approximants to this Borel transform are evaluated, 



V'm(z) 



[M/2] / [(M + l)/2] 



(43) 



surprising that the minimal term in the perturbative ex- 
pansion is reached already in forth order, although the 
coupling assumes the small value g = 0.01. This behav- 
ior immediately follows from the large geometric factor in 
Eq. ( f4l| ) which leads to a "resultative coupling strength 
parameter" of g ros = 0.32. "Nonperturbative effects" of 
the order of exp(— l/g TCS ) provide a fundamental limit 
to the accuracy obtainable by optimal truncation of the 
perturbation series; this is consistent with the numerical 



data in Table III 



We have also investigated the resummation of the di- 
vergent series ( |39| ) via a combination of a conformal 
mapping and Pade aproximants in the conformal vari- 
able. The situation is analogous to the LoSurdo-Stark 
effect: Results are consistent than those presented in Ta- 
ble III obtained by the "pure" Borel-Pade and in this 
case slightly more accurate. The radius of convergence of 
the Borel transform Z' B {z) defined in Eq. p^ ) is s = 1/32 
[cf. Eq. ([I?]) for the LoSurdo-Stark effect], and the ap- 
propriate conformal mapping in this case reads 



4 y 



(i + y) 2 



(45) 



cf. Eq. J22|)|. The inverse reads 



y(w) 



1 - vT=u; 



z' 



[cf. Eq. (Pq)]. The conformal mapping ( |45| ) maps the 
complex w-plane with a cut along (l,oo) unto the unit 
circle in the complex y-planc. While the zero-dimensional 
model example given in Eq. ( |37| ) does not exhibit all 
problematic features of degenerate anharmonic double- 
well oscillators, certain analogies can be established; 
these comprise in particular the need to evaluate the 
mean value of Borel transforms above and below the real 
axis (see Appendix 0). 



where \x\ denotes the largest positive integer smaller 
than x. We then evaluate the (modified) Borel inte- 
gral along the integration contour Co in troduced in jl6| ; 
specifically we define the transform TZuig) 



TZ M (g) 



[ dtexp{-t)P' M 

J Co 



(gt) 



(44) 



In this case, poles above and below the real axis must 
be considered, and the final result involves no imaginary 
part. The particular case of g = 0.01 is considered. 
Values for the partial sums of the perturbation series 
( |39| ) and the transforms defined in Eq. (^ij) are shown in 
Tab. III. The transforms exhibit apparent convergence 
to 6 decimal places in 20th order, whereas the partial 
sums of the perturbation series diverge. Between the sec- 
ond and forth term of the perturbation series, (the forth 
term constitutes the minimal term) , the partial sums pro- 
vide approximations to the exact result. It might seem 



VII. CONCLUSION 

We discuss the resummation of the divergent pertur- 
bation series of the LoSurdo-Stark effect, and of a diver- 
gent model series describing a zero-dimensional theory 
with degenerate minima, using two methods. Method 

I, which uses a variant of the contour-improved Borel- 
Pade technique introduced in pq| , is described in Sec. III. 
The integration contour is modified so that the addi- 
tional terms which have to be added to the principal 
value of the Laplace-Borel integral are clearly identi- 
fied [see also the discussion in J46| and Eq. Jliq)]. Use 
is made of the leading infrared renormalon pole. Method 

II, which comprises an analytic continuation by confor- 
mal mapping with additional improvement by Pade ap- 
proximants in the conformal variable [see Eq. (p~q)], is 
discussed in Sec. IV. This method is a variant of the 
method introduced in [E9lp0| which has been shown to 
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accelerate convergence of perturbative quantum chromo- 
dynamics (by optimal conformal mapping of the Borel 
plane) . We find that both methods accomplish a resum- 
mation of the divergent perturbation series (^) for the 
LoSurdo-Stark effect, and the decay width of the qua- 
sistationary states is obtained (see Sec. |y|. Numerical 
results are given in Tables | and y). A 

main result of the 

current investigation is the demonstration of the anal- 
ogous mathematical structure (doubly-cut Borel plane) 
of the perturbative expansion for the LoSurdo-Stark ef- 
fect and perturbative expansions in quantum chromody- 
namics (renormalon theory |Q ) . The series investigated 
here exhibit a nontrivial singularity structure in the Borel 
plane. In particular, we encounter poles and branch cuts 
on the positive real axis. 

In quantum electrodynamics, we encounter nonper- 
turbative effects in the electron-positron pair-production 
amplitude in a background electric field [72], f^s] [77[ . 
The vacuum-to- vacuum amplitude acquires an imaginary 
part, whose magnitude is related to the production rate 
per space-time volume of fermion-antifcrmion pairs. This 
nonperturbative, imaginary contribution can be inferred 
from the perturbative expansion of the effective action by 
contour-improved resummation (see Jl6[ and the discus- 
sion in Appendix |A|). Nonperturbative effects typically 
involve a nonanalytic factor of exp(— 1/g) where g is an 
appropriate coupling parameter for the physical system 
under investigation (in the case of the LoSurdo-Stark ef- 
fect, the coupling parameter is the electric field strength 
F) . The existence of nonperturbative contributions is in- 
timately linked with the failure of the Carleman criterion 
for a particular perturbation series (see for example [f78[ , 
Theorems XII. 17 and XII. 18 and the definition on p. 43 
in [f79f, p. 410 in |4l|], or the elucidating discussion in 
Rcf. |80[|). The Carleman criterion determines, roughly 
speaking, if nonanalytic contributions exist for a given 
effect which is described by a specified perturbation se- 
ries. 

The current investigation illustrates the utility of re- 
summation methods in those cases where perturbation 
theory breaks down at large coupling. As explained in 
Sec. |v|, even in situations where the perturbation series 
diverges strongly, it can still be used to obtain meaningful 
physical results if it is combined with a suitable resum- 
mation method. In a relatively weak field, it is possible to 
obtain more accurate numerical results by resummation 
than by optimal truncation of the perturbation series (see 
also ) . In a strong field, it is possible to obtain phys- 
ically correct results by resummation even though the 
perturbation series diverges strongly (see the discussion 
in Sec. [v| and the data in Tables | || and [TJ). By re- 
summation, the perturbation series which is inherently a 
weak-coupling expansion can be given a physical inter- 
pretation even in situations where the coupling is large. 
Returning to the analogy to quantum field theory, one 
might be tempted to suggest that physically complete 
results are obtained after regularization, renormalization 
and resummation. 
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APPENDIX A: DIVERGENT PERTURBATION 
SERIES IN QUANTUM FIELD THEORY 

We briefly indicate aspects of certain divergent per- 
turbation series in quantum field theory, in particular 
the quantum electrodynamic (QED) effective action and 
the associated pair-production amplitude for electron- 
positron pairs. We use natural units in which the re- 
duced Planck's constant, the permittivity of the vacuum 
and the speed of light (in ficld-frcc vacuum) assume the 
value of unity (h — eo = c = 1). The one-loop QED effec- 
tive action for an arbitrary electric and magnetic field per 
space-time volume reads [this result can be found e.g. in 
Eq. (4-123) of [[721, upon inclusion of an additional coun- 
terterm; see also [P6lRl|l 



ioo+77 

S = lim -JL f ^ e -(mi-ie )s 




(es) 2 ab coth(eas) cot(e6s) 



-^M 2 (a 2 -6 2 )-l 



(Al) 



The specification of the infinitesimal quantity 77 is neces- 
sary, strictly speaking, in order to avoid the singularities 
of the coth function along the imaginary axis. By a and 
b we denote the secular invariants, 



b= ^T 2 + Q 2 -T 2 , 
T= \f^F^ = ^(B 2 -E 2 ), 

Q = \f^ (*F)^ = -E ■ B . 

If the relativistic invariant Q is positive, then it is pos- 
sible to transform to a Lorentz frame in which E and 
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B are antiparallel. In the case Q < 0, it is possible to 
choose a Lorentz frame in which E and B are parallel. 
Irrespective of the sign of Q we have in the specified frame 

a = \B\ and b = \E\ . 

Then, in the special Lorentz frame, the effective action 
reads [we keep all infinitesimal contributions] 



ioo+7/ 



S= lim -— — 



ds (m 2_ ie);j 
s 3 



(esf \B\ \E\ coth(e|B|s) cot(e|£7|s) 
1 



(esf (B 2 



1 



(A2) 



The particular cases of a pure magnetic and a pure elec- 
tric field follow from the above integral representation by 
considering appropriate limits (\E\ — > and \B\ — > 0, re- 
spectively). These particular cases are of some interest, 
because they can be used as model series for divergent 
alternating and nonalternating asymptotic perturbation 
In the case of a pure magnetic field 



series |]44|-|46 
B = \B\, the result reads [see e.g. Eq. (5) in 



S B 



8tt 2 



ds f , Is 

— { coth s 

s 3 



exp 



eB 



(A3) 



where m c is the mass of the electron. This result can be 
expressed as a divergent asymptotic perturbation series 
in the coupling parameter gs = e 2 B 2 jm\. For the pure 
electric field, the result reads 



of the imaginary part according to the choice of the in- 
tegration contour has been discussed at length in p6|). 
The magnitude of the imaginary part yields the pair- 
production amplitude. The contour C+i is used in the 
current investigation (and in the context of the related 
brief discussion in |4q| ) for the calculation of nonper- 
turbative imaginary effects, i.e. the autoionization decay 
width of atomic states (LoSurdo-Stark effect). 

The divergent asymptotic perturbation series for the 
cases of the magnetic and electric field, generated by ex- 
pansion of the results in Eqs. (A3) and (A4), can be 
found in Eqs. (6) and (7) of [Q (B-field, alternating se- 
ries, coupling parameter gs — e 2 B 2 /m^) and in Eqs. (8) 
and (9) of jl6| (i?-ficld, nonalternating series, coupling 
parameter gE — e 2 E 2 /m^). The singularity structure of 
the Borel transform of the series for the magnetic field 
case has been determined in [ [58) as a sequence of sin- 
gularities corresponding to alternating, factorially diver- 
gent components (these correspond in their mathemati- 
cal structure to the so-called ultraviolet renormalons in 
quantum chromodynamics). The perturbation series for 
the LoSurdo-Stark effect contains both nonalternating 
and alternating components so that its resummation rep- 
resents a comparatively more interesting task. The same 
applies to the more complex perturbation series calcu- 
lated in [|S2| for the renormalization group 7 function, 
whose resummation - at strong coupling - has been dis- 
cussed in f35|,j82f (in this case, there is no imaginary part 
involved). We are not aware of any a priori reasoning 
to determine the absence or presence of imaginary con- 
tributions in a particular physical problem (see also the 
discussion in |83|). 



APPENDIX B: BOREL SUMMABILITY 
IN PROBLEMATIC CASES 
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(A4) 



[Eq. (7) of |4|] and the expression before Eq. (10) of fl|] 
contain typographical errors]. We take the opportunity 
to supplement the pro portionality factor for the expres- 
sion in Eq. (7) of |4(| to yield the effective action per 
space-time volume element; it reads e 2 E 2 /(8n 2 ). As evi- 
dent from the Eq. ( A4 ) , the integration of the Borel-Pade 
transform for the electric field case should be carried out 
along the contour C+i shown here in Fig. [j]. When this 
contour is used, then a sign change results for the imag- 
inary contributions in Table 1 of [[hi (the sign change 



Consider the one-dimensional double- well Hamiltonian 



H{g)=p 2 + x 2 (\~gxf 



(Bl) 



For g — 0, the Hamiltonian describes harmonic oscilla- 
tions. For g > 0, we have degenerate minima of the 
potential at x = and at x = 1/g, and to each eigen- 
value of the unperturbed harmonic oscillator, we have 
to associate two eigenvalues belonging to opposite-parity 
wavefunctions with respect to x = l/(2g). The differ- 
ence of the two eigenvalues is nonanalytic in g. Two 
different approaches have been developed to circumvent 
this problem and to allow for a treatment based on the 
resummation of perturbation series. 

The first approach |6(| is based essentially on general- 
ized Bohr-Sommerfeld quantization formulas and leads 
to an expansion of the ground-state energy eigenvalue in 
terms of 



£( ff ) = E^ (c V + E 



1 



-1/63 



?7=0 
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fe=0 



^fc* 9 ■ 

1=0 



(B2) 



where the upper index of the ^-coefficients labels the 
multi-instanton contributions (the zero-instanton contri- 
bution corresponds to the "usual" perturbative expan- 
sion). The odd-instanton contributions lead to a sep- 
aration of the ground state and the first excited state 
which have opposite parity but the same naive perturba- 



tion series 2~2^Lo 9 l - ^ ne ser i es 2~2^=o ^i"' #' i s non " 
alternating and therefore not Borel-summable; however, 
possible imaginary contributions must be suppressed for 
physical reasons because the ground-state energy is real. 
The suppression can be enforced explicitly by defining 
the sum as the arithmetic mean of the values obtained 
above and below the real axis, or it can be motivated 
by the following observation |60): We define the sum of 



(°) J 



EOO ,-i 
n=0 h l 



C°) J 



g for negative g and carry out an analytic 



continuation to positive g; this leads to an imaginary part 
which in this case cancels with the imaginary part gen- 
erated by the ln(— 2/g) coming from the two-instanton 
contribution (n = 2, k = 1). Note that for the model 



example discussed in Sec. VI, only one alternative is fea- 
sible - the explicit cancellation -, because no additional 
terms are present which could lead to cancellations. 

The second approach |84],[8j| involves contour integra- 
tions and makes use of projection operators in order to 
"select" states of specified parity (this approach has been 
shown to be applicable as well to the problematic Herbst- 
Simon hamiltonian |]86f which involves a vanishing per- 
turbation series). Specifically, we can write the perturbed 
energy eigenvalue as 



E{g) = 



ReFi(g,g) 
ReF (5,sO 



where 



F j (g, 1 ) = -LIdzzH^(g)\ 
2wi J r 



H(rf) 



(B3) 



-I^G?)). (B4) 



where tp^ig) are test functions with a definite parity with 
respect to l/(2g). The closed contour T has radius unity; 
it is chosen to encircle one and only one shifted resonance 
of the perturbed oscillator, while the test functions se- 
lect the state with the desired parity. Specifically, we 
have ip^ig) = -P ± (ff)'0 where tfi is the eigenvector of the 
unperturbed hamiltonian, and the projection operators 
are 



1 



(P ± ( 5 )V) (x) = -W(x)±i,(l/g-x)) 



(B5) 



The functions Fj{g,^j) may be expressed as asymptotic 
series, 



M 

El 

N=0 



j>N + ib j>N (g)} 1 N + 0( 1 M+l ). (B6) 



M+l\ 



The authors of |84, 35| define &f(g, 7) to be the Borel sum 

of X)jv=o a i, N 7 N f° r 9' M an d arg7 small and positive, 
and they establish a corresponding relation for $'((7,7) 
and Y1n=o ^j,N 7 • According to (unnumbered) equa- 
tions on p. 626 of p3, the final result in this case is 
given in terms of the mean values - each obtained above 
and below the real axis - of the two Borel sums and 

' 



ReF,( 5 , 7) = - <f>f(g, 7) + ®f{g,l) 



(B7) 



where z denotes the complex conjugate of z. The value 
KeFj(g, g) is then determined by (unique) analytic con- 
tinuation 7 — * g from RcFj(g, 7). In our simplified model 
example, we have bj.N(g) = (all perturbative coeffi- 
cients are real). The need to evaluate the arithmetic 
mean of Borel sums above and below the real axis appears 
to arise naturally in the context of double- well problems. 
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